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Abstract 

The comprehensive investigation of the temporal evolution of the diocotron instability of the 
plane electron strip on the linear stage of its development is performed. By using the Kelvin's 
method of the shearing modes we elucidate the role of the initial perturbations of the electron 
density, which is connected with problem of the continuous spectrum. The linear non-modal 
evolution process, detected by the solution of the initial value problem, leads towards convergence 
to the phase-locking configuration of the mutually growing normal modes. 
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I. INTRODUCTION 



The diocotron instability is one of the most ubiquitous instabilities Jlj] in low-density non- 
neutral plasmas with shear in the flow velocity. The goal of this study is to gain further 
insight into the physics of the development of the diocotron instability, assessing the role of 
the continuous spectrum 

as 

and the recognition the dynamics of the relative phase of the 
diocotron waves on the linear stage of the instability development. These two problems are of 
paramount importance in understanding the temporal evolution of the initial perturbations 
into patterns peculiar to diocotron instability and has the major practical motivations. The 
solution of both these problems is amenable to examination using an initial value formulation 
of the linearised diocotron instability problem. 

Since the first papers [4| devoted to the theory of the diocotron instability, the com- 
plete analogy of the stability problem of two-dimensional cold low-density electron beams 
in crossed external magnetic and self electric fields, and the problem of stability of the in- 
compressible inviscid shear flows was recognized. Last time, valuable progress was obtained 
in the investigations of the stability of shear flows in a g . A future 

common to these investigations is the recognition that a key element in the understanding 
the dynamics of the instability development, formation of the edges pattern and develop- 
ment of the vortex structures is the relative phase of the edge waves. In this paper we 
adopt this view point to exemain inter-level interaction of edge diocotron waves in what 
is probably the simplest non-trivial case - we consider the plane model of the unbounded 
along the coordinate y non-neutral electron flow, which is confined in the strip —a < x < a 
and moves with ExB velocity in homogeneous magnetic field B = Be z and in own crossed 
electric field, E = E (x) e x , of the electrons strip. Our analysis is grounded on the method of 
the shearing modes (or so-called non-modal approach), developed in first by Lord Kelvin jsj 
(see, also Refs. Jol, [lol]). In the next Section, we formulate the basic non-modal equations. In 
Sec. Ill, these equations are solved for the case, when the initial perturbations of the electron 
density may be ignored and give the modal theory of the diocotron instability. In Sec. IV, 
we consider the diocotron instability in terms of edge waves interactions. The non-modal 
analysis of the diocotron instability, is presented in Sec.V, in which the role of the initial 
perturbations of the electron density on the temporal evolution of the diocotron instability 
to the phase-locking configuration, discovered in Sec. IV, is analysed. A summary of the 



work is given in Conclusions, Section VI. 



II. BASIC EQUATION 



The basic equation in the theory of the diocotron instability in plane geometry is the 
equation for the perturbed electrostatic potential [jj 

Anec dtp drio 



Bq dy dx ^ 
In this paper we consider the homogeneous basic density of electrons, for which velocity 

shearing rate is V^ y = Up e /u ce , where u pe and u ce are electron plasma frequency and electron 
cyclotron frequency, respectively. We consider tenuous electron layer satisfying u pe <C u ce , 
confined in a strip with edges at x — ±a, i.e. (x) = (O (x + a) — © (x — a)), where 
(x — a) is the unit-step Heaviside function (it is equal to zero for x < a and equal to unity 
for x > a); thus, dn (x) /dx = n (6 (x + a) — 5 (x — a)). In this case Eq.([T]) is written in 
the form 



-- Vq(5 (x + a) -5(x-a)) 



dy' 



(2) 



Now we define boundary conditions for potential <p(x,y,t). We suppose, that potential 
is continuous across the surfaces x = ±a, i.e. {x = ±a — e,y,t) = <p(x = ±a + e, y, t) 
with e — > 0. The conditions on the jump of the d<p/dx at x — ±a are determined by the 
integration)^] of Eq.((T]) for the short distances ±e — > across the both surfaces x = ±a, 

d_ ... 8 
dt 
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x=a+e 9X 


x=a—e. 







-Vo-Q-<i>(x = a,y,t) 



(3) 



dt - V °% 



dx 



dx 



x=—a—£-l 



x=—a+e 

-a,y,t) (4) 

Also we require, that potential decays in vacuum region, i.e. {x = ±oo, y,t) = 0. 

We describe two areas: the electron layer, —a ^ x ^ a, and vacuum in the rest of space. 
Eq-© in the vacuum have a form 



dt 



0. 



(5) 



The Fourier transformed over y solutions (x, l,t) — f (ft (x, y, t) exp (—ily) dy of Eq.([5]) are 

<j)(x,l,t) = C x (l,t)e~ lx for x>a, (6) 



(f>{x,l,t) = C 2 {l,t)e 



Ix 



for 



x < —a. 



(7) 



is equal to zero, except the edges at 



In electron layer, the right hand side side of Eq. 
x = ±a, i.e. 

! + ifc£)v, = o. 

In the usually applied normal-mode (modal) approach, the solution to Eq.( 



the form 



(8) 

is sought in 



x) exp (—iout + ik y y), for which Eq.flH]) transforms to the formjj, [U| 

0. (9) 

2, 



(u - kyV^x) V 2 ( 



It is assumed, that u — k y V§x ^ 0, and equation V 2 = is solved|4|, lllj . The proper 
accounting for the resonance u — kyV^x = is attained by solving the initial value problem 
to Eq.([9])j2|, sj- The principal result of the solution of the initial value problem by using the 

was that, in addition to the discrete eigenvalues linked to 
the normal modes, there exists a continuous spectrum of eigenvalues. The application of the 
Laplace transform to the solution of that problem (2, 3] made the calculation complicate and 
low effective - the explicit results for the temporal evolution of the potential were obtained 
only for the asymptotically large time. 

Here we use other approach, which gives easy and transparent treating of the problem 
considered without the application of the spectral transforms in time. It is known, that for 
homogeneous shear flows, the solution of the initial value problem is greatly facilitated by 
a transformation of coordinates x, y to a new set of coordinate frame £,rL that is sheared 
with the mean flow. Such a transformation was proposed by Lord Kelvin [8| and is given by 



(see also Refs. 



X 



£, y = v + v^t, t = t. 



(10) 



In the convected coordinate frame, Eq.Q is spatially homogeneous and the inhomogeneity 
introduced by flow velocity in Eq.([H]) is transformed to a temporal inhomogeneity, 



d_ 

Of 
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d^drj 
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(11) 
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The integration of Eq. (ITT]) over time, 



<9 2 2 y U d 2 (p , ^ , , T ^ 2 \ d : 



oe - / " t d^ + ( 1 + {v " ty JW 

47ren 1 (t = 0,e,ry), (12) 



brings into the considered problem the initial perturbation of the density of electrons in 
layer, n\ (t = 0, £, rj). After Fourier transforming over rj of Eq. (ITS]) , we have 

= 4?reni (* = 0, £, , (13) 

The solution to Eq. (fT3]) for potential <pi straightforwardly gives the initial value problem 
solution for the separate spatial Fourier harmonic for any desired time (without usual 
application [2] of the Laplace transform with time). That solution, 

<P (£, I, t) = (c 3 (/, t) + J dfim (6, e" fcl?1 \ e k ^ 

+ ( C 4 (/, t) + ^J #1*1 (£i, /) e" fe ^ J e fc2 «, (14) 



has an obvious non-modal form with non-separable dependences on time and coordinate in 
exponentials, resulted from the time-dependent k\ y2 = ±1 + iV^lt. The temporal evolution 
of the potential appears to be the strictly non-modal process. 



III. MODAL DIOCOTRON INSTABILITY 



If we suppose that any initial perturbation in layer are absent, i.e. n\ (£i>0 = 0, then 
we obtain the solution, which describes only the surface waves, which form the discrete 
spectrum of perturbations 

0(e,/) = C3(/,t)e^ + C 4 (/,t)e fc2 «. (15) 

The condition of the continuity of the perturbed potential on the boundaries (sc = £ = ±a) 
couples C\ (l,t), C<i (l,t) and C3 (l,t), C4 (l,t), and gives the following presentation for the 
potential: 

(£, = (0> (M) + (I, t) e 2al ) e *- iV ° al \ £ < -a 
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FIG. 1. The dispersion relation of the normal modes. 

(£, = (C 3 (l, t) e'« + C 4 (I, t) e"**) e iV ^ lt , -a < £ < a, 
(£, /) = (C s (/, t) e M + C 4 (/, t) )e-W«« £ > . 



(16) 



Now we apply boundary conditions 
for C 3 (I, t) and C 4 (Z, i), i.e. 

dC 3 



to ffTBT) . and obtain the system of equations 



dt 
dC 3 

dt 



iY±(i-2la)C 3 + i^-e- 2al C A , 
2 v i 6 2 4, 



2 



Vn 

■-±(l-2la)Ci. 



The solution to Eqs. ( ITT)) . (TT8)) gives the relation 



V' 



u = ±-fy/(l-2la) 



2 e -4aZ 



(17) 
(18) 

(19) 



which is the known dispersion equation for diocotron oscillations in plane charge sheet |1 1| . 
The dispersion equation (fl9)) is illustrated in Fig.l. as a function of the parameter al. In 
the case e~ Aal > (1 — 2la) 2 , we have imaginary frequency which defines the growth rate of 
the diocotron instability. In the case, when boundaries of electron layer are so far from each 



other that e < (1 — 2la) , we have two not interacting stable waves. 



IV. MODAL DIOCOTRON INSTABILITY INTERPRETED IN TERMS OF EDGE 
WAVES INTERACTION 

The result of the above analysis is the known dispersion equation (fT9l) . This equation, 
however, does not give an overview of how different initial wave structures will evolve in 
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time. Writing the functions C 3 (l,t) and C4 (l,t) in complex formj^], 

C 3 (l,t)=Q 3 (l,t)e^ t \ (20) 
C 4 (M)=Q 4 (M)e MM) , (21) 

the edge perturbation of the potential can be regarded as two edge waves with amplitudes 
Q 3 (l,t) and Qi(l,t) and phases e 3 (l,t) and e 3 (/,£). By substituting Eqs.( l2"0~|) . (l2"Tj) into 
Eqs. ffT7|) . f[T5|) and separating the real and imaginary parts at x — ±a, we obtain, that 
amplitudes Q3 (l,t) and Q4 (l,t), and the relative phase e = e 3 — €4 of the edge diocotron 
waves evolve according to equations 

^^e-g 3 sm e , (23) 
j t =T (cose + b(t)), (24) 



where 



- + ^ 

and 

, . . 2(1- 2/a) e 2oi 
Q4 Q3 

It follows from (fl9|) . that at the condition of the diocotron instability development, 
(1 — 2la) e 2al < 1. Also, it is easily obtained, that 

(the equality sign appears when Q 3 = Q4). From Eqs. (l2"2"j) . (|2"3"|) one can obtain the integral, 

Ql-Ql = C. (28) 

Therefore, due to the exponential growth of amplitudes Q 3 , Q4 with time from infinitesimal 
beginnings, the amplitudes become almost equal, Q\ = Q\ + C ~ Q\ C, and 6 (t) 
approaches the value &o = (1 ~ 2la)e 2al . Therefore, at the condition, under which the 
diocotron instability develops, we have b (t) < 1, and therefore, the stationary (or fixed) 



points 



12| of the equation fT2"4"|) . where de/dt = 0, exist and are determined by the equation 



cose + bo = 0. The solutions of this equation are two sets of stationary points: stable (or 
attractors) at 

e k = (vr - cos" 1 b ) + 2kn, (29) 

and unstable at 

e k = - (vr - cos -1 6 ) + 2kn. (30) 

Note, that solution of the equation de/dt = cose + bo with initial condition e = eo at 
t — to — 0, has a simple form 



tm £ = -./^f 1 + ^g| (31) 
2 V 1 - 6o I l _ Ae'V 1 ^ ' 



where 



A = (l-6o)tanf + yT^ 
(l-6o)tanf - v^feg 

As it follows from Eq. (l3Tl) . the initial perturbations with arbitrary value of the initial phase 
of each wave, will evolve with time to the ultimate value of relative phase, 



tan l = VT^ (33) 

which does not depend on the initial data and is the same as established by Eg. (1291) (see 
Fig.2). For al = 0.4, Eq.( )33|) gives e ~ 110° (see Fig. 3, where that phase locking configura- 
tion is presented). This solution of the initial value problem reveals the initial stage of the 
instability development as a process of the formation of the phase locked configuration. As 
it follows from Eqs.( l22~]) -(l24"l). the time of the developing of the phase locking configuration 
is comparable with of the inverse growth rate time, t > ~ (Vq) 1 of the diocotron insta- 
bility. Two edge diocotron waves, embedded within spatially distinct regions of oppositely 
directed density gradients interact such, that the wave trains transit to a phase-locked state 
of mutual growth. 



V. NON-MODAL ANALYSIS OF THE DIOCOTRON INSTABILITY. 



In the considered above idealized case, when the initial perturbation rii (t = 0, £, if) of the 
electron density is ignored, Eqs. f|T7|) . f lT8|) resulted from the boundary conditions at £ = ±a, 

8 




FIG. 2. (a)Phase portrait for Eq.(24); (b) The evolution of the relative phase e with time. 

in which non- modal non-separable structure as £i of the spatio-temporal dependence of the 
solution (Il4p receives ordinary modal, as at, form. As a result, we obtain modal theory of 
the diocotron instability. The accounting for the initial perturbation of the electron density 
excludes such simplifications and complete non-modal analysis becomes necessary. In this 
section, we solve for solution (fl~4"]) the boundary value problem, determined by the condition 
of the continuity of the potential 4>i an d by Eqs.Q, (j4j). The continuity of the potential 
gives the following presentation of the solutions (jB]), (0) through the functions C 3 (l,t) and 
C4 (I, t) of the solution (fTSl) in different regions of the space. We obtain for £ < —a 

kia 1 /~i (] ±\ —hia 



(£,M)= [C 3 (l,t)e-^ a + C 4 (l,t)e 

a 

+ / d^ ni (Ci,l)e- kl ^ +a) I e l ^ +a \ (34) 



and for £ > a 



<t>(U,t)= [C 3 (l } t)e k ^ + C 4 (l,t)e k2a 
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FIG. 3. The phase-locking configuration with relative phase e ~ 110°, obtained from Eq, (l30l) for 
al = 0.4. 



27re 



(35) 



The application of the conditions (jlj) to solutions (Il4l . (134|) . and (|35|) gives the inhomo- 
geneous equations for C3 (Z, t) and C 4 (Z, £), 



oc 3 _ .y Q ' 



at 2 



(1 - 2Za) C 3 (Z, t) + i^e~ 2al C 4 (Z, t) + f x (I, t) 



.Vn 



Of 



i^-e~ Aal C 3 (Z, t) - i-£ (1 - 2Za) C 4 (Z, t) + / 2 (Z, t) . 



where 



fi (Z, = / d^m (6, e-^*)* (1 + 2Z - a)) 



(36) 
(37) 



(38) 



a 

/ 2 (Z, t) = -zTre^ I dCmi (£1, eO-*^'*) 6 (1 - 2Z (fc + a)) . 



(39) 



The solution to system fl36l) - fl39l for C 3 (Z,t) and C 4 (l,t) is obtained straightforwardly and 
is given by 



C 3 (l } t) = Cl e^ + c 2 e-^ + C 3 (l,t) 



(40) 



and 



C 4 (Z, t) = ciaie 7 * + c 2 a 2 e- 7 * + C 4 (Z, t) 



(41) 
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where 



4Z 7 

N 2 00 



m=— 00 

x ( e * (-a 2 J! + J 2 ) + <T* (-/ 3 + aiJ 4 )) , (42) 



In Eqs.(J42]), (|43 



Q (Z,t) = £ ni {m,l)j™ 

m=—oo 

x (a^ (-a 2 /i + h) + a 2 e-^ (aj 3 - I 4 )) . (43) 



e -feia=F7ti / 2/ 



= / dti 



a v 1 1 



fciaT7*i / 2Z \ 

- / dh — 1 - 4a/ — , (44) 



gfc 2 aT7*i / 2Z 



'2,4 
t 



- k 2 V - £; 2 



where fci j2 = ±Z + iV 'Zt, 7 = is the growth rate, determined by Eq. (Jl9|) . ai j2 = 
— e 2la ^(1 — 2/a) ± j , and the presentation of the initial density perturbation in a form of 
Fourier series, 

00 

V — "V am? 



was used. Eqs.( TT4|) . (I3"4"j) . (]3"5|) . and ( H2|) -(H51) compose the complete explicit solution on 
the initial and boundary value problems for the separate Fourier harmonic <f>(£,l,t). That 
solution is valid for all times, at which linear theory of the diocotron instability is applicable. 
As it follows from Eqs. (j4"4"]) . ( |4"5]) . at time t < i* = irm/laVl the denominators in integrands 
of and 7 2j 4 decrease with time and approach their minimal value ±Z at time t = £*, which 
is of the order of a few inverse growth rate times. Only after time denominators grow 
with time as t, leading to the decay of the potential with time as (V^'t)~ . It is important 
to note, that only that last stage of thepotential evolution in time is amenable analytically 
by using inverse Lap,ace transform^ g temporal growth of the potentia. at time . < t„ 
known as the Orr mechanism [3, [13J, was completely overlooked. 
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The solution for (£, /, t) for time t > t* is easily obtained by integration of Ji )3 and J 2 ,4 
by parts and may be presented in the form 

0(£,M)=0(o)(£,M) + 0(£,M), (46) 

where 

(o) (e,Z,*) = Ci (e^+^ + a^W) 

+ c 2 (e-^ +fclC + a 2 e-^ t+A;2 «) . (47) 

The constants c\ and c 2 are determined by the initial perturbations of the electrostatic 
potential (or electron density) on the boundary surfaces at x — ±a. The potential is equal 
to 

(f , I, t) = 0! (£, /, f) + 2 (£, Z, + 03 (f , *) , (48) 



where 



0x (e,z,o = c 3 aoe fciC 

/ T/ /\2 oo 

™(V )_ e _ 2la ^ ni(ro>z)e i 



4 7 / 



m=— 00 



fcl" 



nrin 
a 



e kl ^~ u ' I + 



-a) f 



a 2 ai 



\7 + iZaVo' 7 — zZoVq' 

_ e fcl(5+a) n _ 4/fl) f ^ Ol \ 



27 



( 7 2 + (V '/a) 2 ) (Hp - fc 2 ) 
x ( e *i«+«) e -2i« _ (1 _ 4a) e kl ^~ a) e 2la ) } , (49) 



02 (e,Z,t) = C , 4 (/,t)e fc2f 

rra=— 00 



(ssp - Jfe 2 ) V7 - UaV£ 7 + i/a^' 

a 2 ai 
F7 



(inn _ k 2 ) V7 + iia^' 7 - ^a^ ' 
27aia 2 
+ ( 7 2 + (W 2 ) (fc!-^) 
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and 03 is formed by the initial perturbations in Eq.([T 



03 



2vre 



nx(m,l) 



m=—oo 



X 



+ 



i ft- \ iivinE 

e iirm e k 2 {t-a) _ g— ^ 



(Hf - k 2 ) 



(51) 



The presentation of the solution (|46p in coordinates x,y,t is attained by performing the 
inverse Fourier transform for wavenumber I with changing the conjugate variable rj using 
the transformations (fTUl) . 



+0(x,/,t) 



(x,y,t) + 4>{x,y,t) 



(52) 



It is interesting to note, that 0(o) (x, t) has a modal form, 



0(0) (x,y,t) = —J dle Uy { Cl (l)e^ +lx 
+c 2 (I) e^ t+lx + c a (/) a ie ^- te + c 2 (/) a 2 e-^" te ) 



(53) 



whereas (x, t) is non-modal in both sets of variables. The integrand, which determines 
(f>(x,y,t) in Eq. fl52|) . decays with time as (VQ't)^ 1 at times t > and contains non-modal 
multiplier exp (— iV^xtl) (which cancelled in Eq. (l53|) ). From the spectral point of view, 
solution ( |52l) reveals the result, obtained in the asymptotic limit of large time by using 
Laplace transform {2], that the spectrum of the diocotron instability contains discrete 
spectrum of two modes (solution 0( O ) (x,y,t)), and continuous spectrum, which is presented 
by $(x,l,t). 



VI. CONCLUSIONS 



In this paper, we have performed the comprehensive investigation of the temporal evolu- 
tion of the diocotron instability of the plane electron strip on the linear stage of its devel- 
opment. In the realistic laboratory flows, as well as in numerical simulations, perturbation 
growth occurs over a finite time interval due to disruption by instabilities and turbulence. It 
is therefore of interest to find all factors yielding the instability development over a specified 
time interval. We find, that normal mode theory, presented in Sec. Ill, is unable to explain 
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the physics of the evolutionary processes, which occur in the volume and on surfaces of 
the electron strip. The understanding of these inherently non-modal processes attains via 
solution of the initial value problem. By using the Kelvin's method of shearing modes we 
elucidate the role of the initial perturbations of the electron density, which is connected 
with problem of the continuous spectrum 

as. 

Because the growth rate 7 < 0.2VJ, at 
time t > 7 _1 ^ 5 (V^) _1 , effect of the non-modality in the solution for </>(£, l,t) becomes 
important. We find, that linear non-modal evolution process, detected by the solution of 
the initial value problem, leads towards convergence to the phase-locking configuration of 
the growing normal modes, which is the precursor of the development of the multi vortex 
structure, observed in experiments and in numerical simulations. 

It is important to note, that solutions f l42|) -( |45|) . which are valid for any times, incorporate 
additional mechanism of the temporal growth of the initial disturbances, known as the Orr 
mechanism - solutions P2|) - (H5|) experience the temporal growth at t < t*, at which the 
denominators of Eqs. fj4*4"j) . f|45|) achieve their minimal values. Depending on the magnitude 
of the initial perturbations of the electron density, this effect may be comparable with modal 
growth, and it is crucial for the growth of the short along the shear flow perturbations, which 
are stable against the modal diocotron instability. 
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